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The dynamics of a relativistic blast wave propagating through magnetized medium is considered 
taking into account possible inhomogeneities of density and magnetic field and additional energy 
supply. Under the simplifying assumption of a spherically symmetric explosion in a medium with 
toroidal magnetic field self-similar solutions for the internal dynamics of the flow are derived. In 
the weakly magnetized case, when the bulk of the flow may be described by the unmagnetized 
solutions, there is a strongly magnetized sheath near the contact discontinuity (when it exists). Self- 
similar solutions inside the sheath are investigated. In the opposite limit of strongly magnetized 
upstream plasma new analytical self-similar solutions are found. Possible application to the physics 
of Gamma-Ray Bursts is discussed. 



1. Introduction 

Recent observational advances in the area of Gamma- Ray Bursts (GRBs) (see, e.g., Ref |^ for a review) have stimulated 
research on the physics of strongly relativistic explosions. Most analytical results for the dynamics of strong explosions 
belong to the class of self-similar solutions. This approach allows great simplification in reducing a system of partial 
differential equations ia the ordinary ones and often represents an asymptotic behavior af the flow. Following the 
seminal works of Sedovo (for non-relativistic strong explosions) and of Blandford & McKeed (B&M hereafter) a series 
of generaliziug works have been done. Non-relativistic shock waves propagating in a magnetized medium have been 
consideredoO as well as generalizations of the B&M solutions for a wide class of density profilcsLrU. Yet, until now 
no account of the dynamics of the magnetic field in relativistic blast waves has been made. A lack of the relativistic 
treatment of the blast wave in magnetized medium is noticeable, especially in view of heuristic interpretation of 
GRB emission as due to synchrotron emissiouJn strongly magnetized relativistic blast waves (e.g., Ref This is 
done in the present work. In a follow-up workllJ we will explore the internal structure of the relativisticly expanding 
magnetic cavity. 

2. Formulation of the problem 

We seek self-similar solutions to the relativistic dynamics of a spherical expansion of gas into a magnetized medium 
due to strong explosion. We assume that (i) time dependent energy source is located in the center; (ii) the magnetic 
field is toroidal and spherically symmetric; (ii) there is a cold spherically symmetric external medium with density 
ni (which may depend on radius). 

The assumption of spherical symmetry with toroidal magnetic field is a controversial, but |sUGommonly employed 
simplification. Ajrmmber of works applied spherical approximation to MHD winds from pulsarscJilj and to explosions 
in the Solar windQ'Q. Other astrophysical setting where such situation may arise is an explosion in a preexisting cavity 
blown out by spherical magnetized stellar wind. Formally, our approach also applies to the equatorial region of an 
explosion in a constant magnetic field, where the shock velocity is perpendicular to the magnetic field. 

The formal treatment of the problem, in which we rely mainly on the approach of B&M, starts with the set of 
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F;y=0, (2) 

{pu')..^0 (3) 
where 

T'^ = {w + h^)u'u^ + {p+ \b'^)g"' - b' V (4) 

is the energy-momentum tensor, w = 4p is plasma proper enthalpy (upfront in the cold plasma wi = nimiC^, where 
ni is external density, rui is the ion mass, while behind the shock w = 4p appropriate for relativistic plasma), 
6^ ~ biV is the plasma proper magnetic energy density times 47r, p is pressure, — (7, 7/3) are the plasma four- 
velocity, Lorentz-factor and three-velocity, g^^ is the metric tensor, 6^ = ^rjiji^iu^ F^^ are the four-vector of magnetic 
field, Levy-Chevita tensor and electro-magnetic field tensor. It has also been implicitly assumed in the derivation of 
these equations that one of the electromagnetic invariants is not equal to and the electro-magnetic stress energy 
tensor can be diagonalized (equivalently, this implies that there is a reference frame where the magnetic (or electric) 
field is not equal to 0). 

Writing out Eqns in coordinate form and assuming a spherically symmetric outfiow with toroidal magnetic 
field, the conservation of energy and momentum (^, induction equation (|2|) and mass conservation give 

dt [{w + b^)j^ -ip + 52/2)] + ^dr [r^{w + b^)f3-/^] = (5) 
dt [{w + 6^)72/5] + [r^ ((^ + 62)^2^2 _^ ^ ^2/2))] - ^ = (6) 

dt [67] + -dr [rbp^] = (7) 
r 

dt [P7] + \dr [r^pPj] = (8) 



Equation (g-j^) should be complemented by the boundary conditions on the shock front. The jump conditions on 
the relativistic transverse magnetized shocks may be written aal3 

_-r2 ( h if cr < 1 



7 = 


7' + u' 


P2 = 


m ni 
4u'7' 


n2 = 


niT 


62 = 


u' 


where 





'il 

:..(:- 22;)) 

J 2V2 if < 1 

r / 9 / f 2 \/2ct if (t < 1 

1 if cr > 1 



niF X 



(9) 



u 



,2 1 + 10(7 + 8^2 ^l + 20a (l + a) + 64a2(l + a)' ^2 ,2 



16(l + cr) 16(1 + ct) 

are the four speed and Lorentz factor of the shocked fluid in the frame of the shock and 
J, 2 
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is the magnetization parameter in front of the shock. The subscripts denote the quantities in the unshocked (1) and 
shocked (2) media. 

Relations (p|) allow for the following parameterization 
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P2 = V^nimc^ f{x) 
n2 rnin(x) 

62 = Th^Kx) (12) 

with the boundary conditions /(I) = /o, 5(1) — go,h{l) — ho, n(l) = no (see below). ^ 
Following B&M we choose the self-similar variable 

X = 1 + 2(m + 1)^ = [1 + 2(m + IjT^] (1 - r/t) (13) 

where = (1 — r/R)T^, T is the Lorentz factor of the shock, R = t{\ — l/{2{m + l)r^)) is the radius of the contact 
discontinuity and we assume the Lorentz factor scales with radius as cx t^™. We limit ourselves to the strongly 
relativistic case expanding all relations to first order in 1/F^. 

Treating (Xi ?/)j where y = F^, as new independent variables we find 

dt = -mydy + ((m + \){2y - x) + l)^^ 
dr = -{l + 2{m+\)y)d^ 

-'('' 1 + 2(^1), ) ("> 

The equations for the self-similar variables /, g, n and h read 

Adlnf , ^ , 2(l-m) + 3 (m-2) xg ,2 

^ — = -2 4 1 - m + m - 4) X .g) / + r tt 
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^ 91n/i _ (2 (to-4) x^9^ + (8-llH (xg-1)) 3 (to - 1) ^^2 
g dx ~ 2(x<?-l) -^+2 

= TT (-.^ (-12 + 11 ni + g (24 - 11 to) x + 2 5' (to - 6) x') + 

9 ox 2(x.g - 1) ' ' 

ih? (1 — TO + g (m — 3) x)) 

^=(i + to) (2/ (i-4x.g + g'x')-3x5^') (15) 

System ( p^ ) is a relativistic generalization of Eqns. (3.10-3.13) of Ref. ^ In the limit h ^ {) these relations reproduce 
the unmagnetized case of B&M. 

Sometimes it is more convenient to use use the functions g,f,h,h and x — gx with boundary conditions on the 
shock front g~f = h = h = x = l: 

9 = 9q9, f = fof: n = noh, h = noh, x = gx ^ gox 



Y -U' [j^ cr > 1 



^We choose to work consistently with proper quantities, i.e. measured in the plasma rest frame. One should be careful in 
comparing our equation with B&M and Ref. ^ 

^ The system (^ may in fact be transformed into a system of one ODE and three equations resolvable in quadratures if we 
change to a new coordinate x — xd and introduce a magnetization parameter /3 = f/h^- 
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We can solve the above equations by direct numerical integration, (see Fig. [ijand Simple analytical relations, 
though, can be obtained in the weakly, cr <C 1, and strongly magnetized, cr ^ 1, limits. 

3. Small a limit 

Relations for g,f and n in this case are the same as for the unmagnetized case (B&M). For the evolution of the 
overall passive magnetic field we find 

Idlnh _ (-16 + 22TO+ (8- 11m) X5+ (m-4) x^ff^) 

^~d^^ 2 (1 + m) (x5-2) {i-8x9 + X^P) ^ ^ 

which can be integrated using the known solution for g. In particular, for m = 3, this can be resolved for h{x)' 

Kx) = I (18) 

A 

Recall that in this case g = 1/x, f = 1/x^^^^^ and n = Xjx'^^- The magnetization parameter, /3 = j 1^ = 1/x^^^^ 
in this case is a decreasing function of X; so if it is small on the shock it will always remain small. 

A. Magnetosheath 

When the energy contained inside the shock increases with time a power source is necessary. This can be an expanding 
piston or another fluid with a larger Lorentz factor. In both cases a contact discontinuity (CD) forms separating the 
external and internal fluids. As discussed by B&M the point x<? = = 1 is the location of the contact discontinuity. 
As Eq. (0) (which is the small a limit of exact relations) suggests, the magnetic field increases infinitely near the 
CD. In reality, of course, it will grow until approxirriaijC equipartition is reached forming a thin boundary layer. This 
is consistent with the non-relativistic considerationll£lEl where it was noted that however small the magnetic field is 
in the unshocked medium near the CD the magnetic field grows in magnitude and starts to dominate the dynamics 
of the thin boundary layer. In this section we study the dynamics of the magnetosheath. 



To estimate the thickness of the magnetosheath we use the system (A2) using independent variable a; (fc = / = 
in this section). In the small cr limit the local magnetization parameter /? — / f is 

fi ^ {X- l)2(™-4)/(12-r«) (jg) 

Since near the CD g ~ constant the width of the magnetic sheath in coordinates x is 

Ax OC a(12-™)/2(4-m) (20) 

The constant of proportionality here is a complicated function of m. The width of the magnetic sheath in coordinates 
X is 

Ax = , (21) 

9cd/9cd + gcD 

where g'fjjy and gcD are the values taken at the CD. 
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Next we study the structure of the boundary layer. Using x as independent variable and retaining leading terms 
near a; = 1 we find 

dlnf (to -4) 



dx (6/i2 _ / (m- 12)) {x - 1) 

d\ng _ 3/i2 (to-1) + / (-8 + 5 to) 



dx 


-6/i2 + / (to - 12) 


d\nh 


/ (to - 4) 


dx 


[Gh^-f (to -12)) (x-1) 
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(2/to + 3/i2 (-1 + 3x)) 



(22) 



dx 2 (6/i2 - / (to - 12)) {x- 1) 
which immediately gives 

/ + hy2 = Ci (23) 

where Ci determines the total energy flux through the magnetosheath. For simplicity below we put Ci equal to 1. 
Using ( ^3|) we can integrate the equation for the magnetic field 

^™-i2 ^^2 _2f{x- 1)'"-^ = C2 (24) 

where C2 is a constant of integration. Relation (p^ ) implicitly determines the self-similar structure of the magnetic 
field inside the boundary layer. 
As we approach the CD, a:; — *■ 1, 

/i - V2 - (a; - l)(4-")/6 (25) 

This shows that for to < 4, as a; ^ 1 the magnetic field goes to a constant on the CD: h \/2. Similarly we find 
that f ^ {x — l)(''~™)/6 and n ~ y^x — 1 0. Thus, pressure and density vanish at the CD while the magnetic 
field is finite. This result confirms that no matter how small the external magnetic field is, there will be a region 
near the CD where magnetic field pressure is dominant. This is qualitatively different from the hydromagnetic case 
where (for m > 0) the density was vanishing on the CD, resulting in very high temperatures (c.f. Ref 

The assumption of ideal MHD should break down somewhere inside the boundary layer where diffusive and 
dissipative processes are likely to play a very important role. 

The solution inside the boundary layer should be matched to the solution in the bulk flow. The exact relations 
connecting Ci and C2 to a and to are prohibitively complicated for reproduction here. 



4. Large a limit 

Simple analytical results may be obtained in the case of a strong magnetization, 1 ^ a ^ V'^ /A. The upper limit 
on (T comes from the fact that for a >r^ /A the shock is no longer strong: it's four-speed become comparable to the 
upstream Alfven four-velocity. 

Changing to tilde-functions and expanding Eq. (p^ ) for cr S> 1 we find 
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1 91n/ ^ 2h^ (m-1) ^^e) 
9 (1 + m) (/-3/i2xg) 

Which immediately gives 

9 = h^, f = h'^'^, n = h (27) 

In the strongly magnetized limit the pressure / is proportional to density n to the 4/3 power. 
We can then resolve the equations for h as an implicit function of x- 

2 (3 — m) - 2 (i+,Ti) 1 -|_ TO g 

X = ^ r^h-^^ + -h~^ 28) 

(7 — m) (7 — m) 

(see Fig In particular, for the point explosion case, m = 3, this gives 

h = x-'^\ ~f-X-"\ ~9 = X~"\ n^x-'" (29) 

The local magnetization parameter [3 — h? / f = X~^^'^ is a slowly decreasing function of x- 

Relations (E9|) represent new self-similar solutions for strong point explosion in a strongly magnetized medium 
with 1 < (T < r2/4. 

As a test to these solutions we note (following the arguments of B&M) that for point explosion the energy associated 
with some interval of dx should remain constant. This leads to the condition 

where (3n = (ff)^-consf leading orders of 1 <C cr ^ this requires 



which is satisfied by the relations (]2£ 

The total energy contained inside the relativistic strongly magnetized shock is, for m — 3, 

r" TT 

E w int^ / j'^h'^r'^dr -t^nirriiC^y (32) 
Jo 4 

which is independent of time. This relation provides the normalization for y — T^. 
For the case of a blast wave with an energy supply we find: 

E = C{m) nim,c3i3y4/(i+m) ^ i(3^m)/(m+i) ^33) 

with the constant of proportionality C{m) being a complicated hmction of to. If the energy source is a power law 
function of time L — L^f^ , then 

m = 1^ (34) 

(c.f. B&M Eq. (57)). 



5. Blast wave in an inhomogeneous medium 

In this section we assume the unshocked density and magnetic field have power law dependences on the radius: 
ni ~ and hi ^ r^K As a classical example the non-relativistic constant velocity wind gives k = 2, I = 1. 
Straightforward calculations give 
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= (4 + m ixg - 4) + fc (X5 - 2) - 4x<?) + 
ft.^ (-2 (/ + m - 1) + (-6 + 3 fc - 2 Z + 3 m) X5) 

/ (-4 + 3fc + 7m-2 (2 + m) X5) + 3 /i^ (l + m - 1) 

9 ox 

Adlnh _ f {-8 + 3k + 4:l + Um+{8 + 3k-16l-nm) xg + '2 {-4: + 2l + m) x^9^ 
9 dx 2 {xg - 1) 

3h'^ jl + m-l) 
2 

^Slnn _ / (12- 7fc- llm+ (-24+ 13A:+ IIto) x.9 - 2 {-6 + 2k + m) x^9^) 
9 dx 2 (x5 - 1) 

3 /i^ (/ + m - 1 + (3 - 2 /j + / - to) X5) 



2 (X5-1) 

^=(l + m) (2/ (l-4x.g + x'.9')-3x5/i') (35) 

In the limit of small a we reproduce equations (62-64) of B&M {I naturally falls out) plus an equation for the 
magnetic field: 



(36) 



Idliih _ (-16 + 6fc + 8; + 22m + (8 + 3fc- 16 Z- 11 to) x5 + (-4 + 2Z + m) x^ 9^) 

2 (1 + m) (xg-2) (4-8x5 + x'5') 
which for the point explosion case m = S — k (magnetic energy density is not important) gives 

^^i, h = x''^^, /^;^(17-4fc)/(3/c-12)^ -^^ (10-3fc)/(4(/c-4)) (37) 

X 

If the upstream medium has a constant magnetization parameter, then I = k/2 and h — 1/x- 
In the large a limit we find g = h? , / = ft,^/'^, h — h and an equation for the magnetic field 



Idlnh 3h{l + m-l) 



9 9x 2 (1 + to) (/-3/i2 gx) 
which may be resolved for x(^) • 

1 / , - 2(ll-8i-5r,i) 



(38) 



X 



7 ~2l-m 



/ ^8 , 2 (11-8 i-5 m) \ 

((1 + to);i"3 +2 (3-2Z- m)/i 3(i+™-i) j (39) 



Again, similarly to the homogeneous case, the point explosion, m — 3 — 21 corresponds to h = 1/x'^^* The case of 
constant energy source in a current-free wind (I = 1, m = 0) gives g — f = h — n — I 

Similarly to the homogeneous case we can obtain the solutions inside the magnetosheath for the case of nonzero 
k and I. We find that the width of the boundary layer in coordinates x is 



12-3 k- 



AX = (7 2(-4 + 3fc-4 1 + ™) (40) 

The structure of the magnetic field in the layer is determined implicitly as 

a: - 1 = C2 (2 - /t^) "-■+^'=-'^'+"/t 4-wit". (41) 

(compare with (|^). Relation (^) determines the behavior of the magnetic field on the contact discontinuity. For 
example, for a constant external magnetization parameter I = k/2 magnetic field completely dominates the boundary 
layer dynamics for fc < 4 and fc + to < 4. For a shock wave propagating in a non-relativistic constant velocity and 
constant magnetization wind (fc = 2, Z = 1) magnetic field goes to a constant while density and pressure go to zero 
for to < 2. 
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6. Discussion 



We have analyzed the propagation of a spherically symmetric relativistic shock wave into an inhomogeneous medium 
permeated by toroidal magnetic field. The main limitation of the current work is the neglect of the possible lateral 
expansion. The natural and much more complicated extension of this approach would include a blast wave in a 
constant magnetic field. 

The self-similar solutions presented here are generally applicable to the expansion of a strongly relativisticly blast 
wave in a preexisting spherical cavity blown by the magnetized wind and in a equatorial region of a constant external 
magnetic field. The self-similar structure of a thin magnetosheath, on the other hand, has a much broader validity 
since any radial (orthogonal to the surface of the contact discontinuity) magnetic field should approach on the 
surface of the CD. 

The results of this work bear relevance to the physics of Gamma Ray Bursts. In particular, the most ijMjpular model 
of GRBs relies on relativisticly strong shock waves to produce both Gamma-ray and afterglow emissioncl. Previously, 
the dynamics of the magnetic field has not been included in the model. In addition, the very basic assumptions of 
the model - prjesence of near-equipartition magnetic fields and acceleration of particles by relativistic shocks - have 
been criticizedEJ. Our results point to the new interesting possibility that the relativistic flow may produce radiation 
effectively in the magnetosheath layer adjacent to the contact discontinuity between two flows where magnetic field 
reaches near-equipartition values. We leave the investigation of this possibility to a future work. 
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Appendix A: Dynamic equations for independent variable x = gx 

Changing to the new variable x = gx in ([l5| ) and using 

{-iK^ {-2 + 1) x + f {-2 + m{x-2)-i{k-A) x + 2x'^)) g{x) 
'^'^ " (1-fm) {ih^x-2f {l-^x + x-')) '^^ 

we find 

(x-1) A! 2f U + m{-A + x) + k {x-2)-Ax) ix- 1) 

ox 

(2 (/ + m - 1) + (6 - 3 fc + 2 i - 3 to) x) 

A! f U--ik-7m + 2{2 + m) x)--ih^ n + m-l) 

ox 

ix-l)M^^^^ -f{-S + ?,k + Al + nm+{^ + ik-\Ql-llm)x + 2(-A + 2l + m)x^)- 
Ox 2 

hi^ [l + m-l) [x-l) 

d In Tt 1 

ix - I) A! — — = - (-12 + 7fc + llm+ (24- 13fc - 11m) a; + 2 (-6 + 2fc + m) a;2) / + 
ox 2 ^ ' 
3 

- (? + TO - 1 + (3 - 2 A; + ; - to) x) /i^ 

-3 (-2 2; + / (-2 -I- TO (x - 2) - 3 (fc - 4) x -I- 2 x^) (A2) 
In coordinate x the shock is located at xq — go and the contact discontinuity is located at x = 1. 
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Appendix B: Lower dimension outflows 



Without giving detail calculations, we comment here on the generalization of the above results for the lower di- 
mensional systems, when the outflow is cylindrically or linearly symmetric. We give here only the scalings for the 
point explosion cases in a homogeneous medium, generalizations for inhomogeneous media with energy supply are 
straightforward. For 1-D case (strong shock propagating along a tube) the solutions look very similar to the 3-D 
case. In the small a limit we find g = h ^ n = 1/x, / = X^^^^, while for the large a limit we reproduce relations 
(p9|). The structure of the equations in the 2-D case is, in fact, qualitatively different from the 1-D and 3-D cases. 
In the small a limit: g = h = l/x^^^, = 1/x^^^i / = X~^^^j while the large a turns out to be unusual: 

~g = ;t-5/6^ fi ^ h ^ X^'^'', / = X'^'^- 
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Fig. 2. Lorentz factors, pressure, magnetic field and density for m = 1. The curves are labeled by the values of a. For 
cr » 1 the asymptotic solutions are g = f = h = n = l. Note that magnetic field piles up on the contact discontinuity. 
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